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■ ABSTRACT 



We investigate the electromagnetic fields in the vacuum exterior of a rotating relativis- 
tic star endowed with a magnetic dipole moment, and with the stellar surface behaving 
as a perfect conductor. While the stellar rotation is treated in the slow-approximation 
of general relativity, we do not restrict our attention to slowly rotating electromagnetic 
^) , fields, and take our analysis beyond the low- frequency approximation considered so 

£SJ ' far. When the dipole moment is misaligned with the rotational axis, our approach does 

not yield analytic solutions as in Rezzolla et al. (2001a, b), but determines the prop- 
erties of the electromagnetic fields approximately and semi-analytically, by computing 
the coefficients of simple expressions for the fields through the numerical solution of 
two partial differential equations. Because our approach provides a solution which is 
in principle valid throughout space, we can evaluate the accuracy and/or invalidity of 



Q ' previously known analytic expressions at different distances from the stellar surface. 

, Overall, the solutions found in this way represent an efficient way of bridging in a 

^ ■ single semi-analytic formalism the strongly relativistic (Rezzolla et al. 2001a, b) and 

the asymptotic regimes (Deutsch 1955) for which analytic solutions have been found. 
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1 INTRODUCTION 

It is well known that magnetic fields o ften play a key role in astrophysics and an ample discussion is given, for instance, 
in IParke"rl l|l97flF ): Tzeldovich et alJ l)l983|) . In the limit of very high electrical conductivity, the magnetic flux is "frozen" into 
the fluid and increases as a result of fluid compressions. In this case the the typical value of the magnetic field B can be 
related to the rest- mass density p through the simple relation B oc p 2//3 , so that very strong magnetic field may be produced 
in compact stars, where the gravity is also strong. Some observations indicate that the magnetic field in young neutron 
stars is of order 10 12 G. Recently, exceptionally str onger cases with 10 14 - 10 15 G are suggested in the ph enomena with 



soft gamma repeaters and anomalous X-ray pulsars ijMereghe tti g^U^||l99jj; [Kouye liotou et al1ll998l . Il999h . Very intense 
magnetic fields may be produced through a dynamo action (Tho mpson fc DuncarJll993f) and the neutron stars with such 
magnetic fields are usually referred to as magnetars. The ideas behind the existence and formation of magnetars are gaining 
support from the observations and a number of authors have discussed the mecha nisms generating these magnetic fields and 
their astrophysical relevance (see, for instance. iThompson fc Duncan! l)l99.4 Eifflrt lV The combined effects of strong gravity 
and intense electromagnetic fields may be crucial in some astrophysical objects and it is therefore important to consider 
the electromagnetic field in curved spacetime produced by strong gravity. The study is important in clarifying the general 
relativistic effects on the electromagnetic fields, and also indispensable in constructing realistic models. 

A number of works have been so far published concerning the effects of strong gravity on the electromagnetic fields. The 
magnetic field in principle affects the spacetime curvatu re through the Einstein e quations, but the actual deformation at the 
surface is of order l(T 5 (B/lu 15 G) 2 even for magnetars (|Konno et al.l In most astrophysical cases, therefore, the 
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deformation may be neglected, and general relativistic effects can be evaluated by solving the Maxwell equations in fixed but 
curved spacetime. Several attempts are done along this line of t hough t . Stationary electrom agnetic fields have been considered 



curvea spacetime. several attempts are done aiong trus line or tnougnt. stationary electromagnetic ncias nave dccii consiaerca 
both in Schwarzschild spacetim es (e.g. iGinzburg fc Qz ernoy I (Tl964h : lAnderson fc Cohen"! 1 1970|}JPe^tersonl (h97jh) and in 
slowly rotating spa cetimes (e.g. iMuslimw' fc'^Vs ganTilTflOSl) : iMiislhrroT ^TlfardiT^T^lQQ? ! : KotmolTl^olmial jgflflfl) ). 

More recently, iRezzolla et al.l l|20^^aH bfl and lzaiiottifcRfizzoll^l(|200^) have extensively studied the interior and and 



vacuum exterior electromagnetic fields produced by a rotating dipole moment comoving with the star. Their analysis is fully 
general relativistic and has provided compact and analytic expressions for both the cases in which the magnetic dipole is 
aligned or not with the stellar rotational axis. In these works, the star of mass M and radius R is assumed to be rotating 
at an angular velocity f2* much smaller than the maximum angular velocity, i.e. the Keplerian angular Q K ~ -J GM/R 3 , so 
that a slow-rotation approximation can be used satisfactorily to describe the general relativistic corrections introduced by the 
rotation of the spacetime. Furthermore, the electromagnetic fields are assumed to be comoving with the star at all times, thus 
possessing an intrinsic angular velocity £7 = Q, -C Q K ; as a result, the expressions of Rezzolla et al. (2001a, b) are expected 
to be accurate in the vicinity of the star, where the general relativistic corrections are stronger and the intrinsic velocities 
smaller. We will refer to this as to the low-frequency approximation in the Maxwell equations. 

While the slow-rotation approximation is a good one except for extremely rapidly rotating stars, the low-frequency 
approximation may have some limitations. The first one is that the expressions of the electromagnetic fields derived in this 
case cannot be accurate at larger radii, (i.e. in the vicinity and beyond the light cylinder), where they should manifest a 
wave-like nature. The second limitation may arise in the case in which the electromagnetic fields are not comoving with the 
star and, more specifically, in the case in which Q ^> fi„. 

In this paper we circumvent these limitations by considering the general relativistic electromagnetic fields exterior to a 
rotating neutron star in the slow-rotation approximation but not in the low-frequency approximation. We will not discuss 
here the physical conditions that may lead to electromagnetic fields with Q 3> fi*, nor whether this is likely to happen 
for realistic astrophysical sources (indeed this condition may be difficult to occur in practice); rather we will make this our 
working hypothesis. As a result, we here consider the electromagnetic fields in the vacuum exterior of a relativistic star 
slowly rotating at frequency f2* <C fl K and endowed with a magnetic dipole moment rotating at a frequency (either Q = f2* 
or fi ^> fl« ). Furthermore, the stellar su rface is assumed to be a perfect conductor, thus inducing a quadrupolar electric 
field jDeutschlll95llRezzolla et alJbOQlal) . Our approach differs from previous ones both in the underlying assumptions (i.e. 
high frequency electromagnetic fields) but also in the method. In particular, we do not provide analytic expressions (as in 
Rezzolla et al. 2001a, b) but determine the properties of the electromagnetic fields through a semi-analytical approach in 
which simple expressions for the electromagnetic fields need to be completed through the calculation of coefficients coming 
from the numerical solution of two partial differential equations. Because the approach is based on a series expansion, the 
solution and its accuracy depends on the number of terms considered. On the other hand, our approach yields solutions that 
are valid in the whole spatial domain and this allows us to compare our results with the analytic expressions of Rezzolla et 
al. (2001a, b) and of Deutsch (1955), estimating the accuracy of these analytic solutions in those regimes where they are not 
expected to be valid. 

The paper is organized as follows. In Section 2 we summarize the Maxwell equations in vacuum with stationary condition. 
The equations are reduced to a set of two partial differential equations among two scalar potentials. Here, we discuss how 
it is possible to to look for simple Deutsch-type solution which are approximate solution of the general relativistic Maxwell 
equations in a slowly rotating spacetime. In Section 3, we show the explicit calculations in the case in which the magnetic 
dipole is parallel to the rotational axis, while in Section 4, we consider the case when the dipole is perpendicular. Because 
of the linearity of the Maxwell equations, the general case of an oblique rotator can be constructed from the superposition 
of these two particular cases. Concluding remarks are finally presented in Section 4. Hereafter we will use units in which 
c=G = l. 



2 STATIONARY ELECTROMAGNETIC FIELD 

We will consider the solution to the Maxwell equations in the fixed background spacetime of a rotating relativistic star whose 
generic axially symmetric line element is given by 

ds 2 = -adt 2 + jijidx* + Pldt){dx 3 + Pldt) (1) 
= -a 2 dt 2 + e 2 ^ {d<t> - Ludt) 2 + e 2 "dr 2 +e 2 "d(9 2 . (2) 

In practice and except for rapidly rotating stars, it may be sufficient to take into account the first-order rotational effect of 
the metric only (see, for instance jZanotti fc R,ezzollall2002T) 1 so that the metric functions exterior to the slowly rotating star 
are reduced to the form 

(, 2M\ 1 ' 2 „ / 2M\- 1 /2 . 2 J 

a — v 1 , e — ► 1 , e — ► r, e — ► r sm a, u> — > — =- 



r J \ r ) 



r-3 



(3) 
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where M and J are the mass and angular momentum of the central star. In this paper, the electroma gnetic fields are denoted 
by the values measured by "zero angular momentum observers", i.e. ZAMOs feardeen et alJll972ft . and their components 
projected to a locally orthonormal tetrad frame are expressed as (E r , Eg, E^) and (B r , Bg, B$) . These are sometimes expressed 
with hatted indices to distinguish them from coordinate values (see, for instance, Rezzolla et al. 2001a, b); however, because 
we will always refer only to the electromagnetic fields as measured by ZAMOs, we will omit here the hatted in dices. 

The electric and magnetic fields in vacuum satisfy the Maxwell equations in this curved spacetime (e.g. iThorne et alJ 

<198dn . 

V • B = 0, (4) 
V-E = 0, (5) 
d t B - L^B = -V X (olE), (6) 
d t E - Lp E = V x (qB), (7) 

where Cs is the Lie derivative along the vector field /3* = —uidj, and is defined for any vector A as 

C 0f A = 0*-V)A-(A-V)p*. (8) 

Clearly, in the frame rotating with angular velocity j3 = Q.d^, the electromagnetic fields are independent of time and these 
are the electromagnetic fields we will consider hereafter. Note that the angular velocity fl may differ from that of the star 
fl„, which is related to the spin angular momentum J and the inertial moment / through the relation fi„ = J/I in the slowly 
rotating spacetime. For example, the electromagnetic fields can be produced by a periodic motion with a certain angular 
frequency Q on or outside the star. Furthermore, and as mentioned in the Introduction, we do not assume that the frequency 
O is small, unlike in previous work. The stationarity of a scalar function / is equivalent to expressing that / should depend 
on the time t and on the azimuthal angle <f> only in a combination of the type </}' = (/) — fit. In the locally flat four-dimensional 
spacetime in which the ZAMOs make their measurements, this condition of stationarity for the function / can be expressed 
by requiring that the function / is solution of the equation: dtf + 13 ■ V/ = 0, i.e. that the function is a constant along the 
direction (l,/3), in the four-dimensional reference frame of the ZAMOs. In the case of the electromagnetic vector fields E and 
B, the mathematical condition for stationarity needs to be expressed through the Lie derivative as 

d t E + C E = 0, d t B + C B = 0. (9) 

Using these conditions in eqs. © and Q, some combinations of E and B are expressed by gradient of two scalar potentials, 
H and G. After solving for E and B, we have 

E= — V le-^Cr +^-"a~ 1 J/, e , e- v G,e -^'V 1 ^ , e^AGU 1 , (10) 

B — — V \e~"H,r -zue^a^Cg , e' v H, e -f^e^ci^G,,. , e^AH,A 1 , (11) 
qA l j 

where a comma denotes a partial derivative with respect to the coordinate j(= r,0,<j>) and 

A = 1 - ^e^a- 2 , (12) 
zu — Q — uj. (13) 
From eqs.0 and ©, we have a pair of second-order partial differential equations 

A a (e- f ** + *ar 1 G,r) r + a (ef-*+*ar l G,o) g + e" +, "*AG^^ 

^+i>^ n ^~^A,gG,g + K r H,g-K e H, r = 0, (14) 

,o) g + e^-^AH^ 

A, r H, r - e"~ v+ ' p A, e H : g - K r G,g + K e G, r = 0, (15) 



A 



where 



M . (for j = r,6). (16) 



,3 



These equations 1141 and 1151 are symmetric under the transformation between E and B, i.e. E — > B, by changing G — » H 
and H — > — G. Note that this is a generic property of the source-free Maxwell equations. 
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AXIALLY SYMMETRIC CASE 



We now look for solutions of eas. 1141 and 1151 in the slowly rotating spacetime. Note that these equations have a singularity 
at the light cylinder A = 0, so that a special care should be paid when performing a numerical integration. We now expand 
the electromagnetic fields using a series expansion with respect to the angular part which is particularly advantageous for 
the imposition of regularity conditions. Since the electric and magnetic fields are respectively polar and axial vectors with 
respect to parity transformation, we expand them through two scalar potentials G and H, which have parity of +1 and — 1, 
respectively 1 . Therefore, for an axially symmetric system, these functions are respectively expanded in terms of the spherical 
harmonics an d Y^+i- Note that the following forms are equivalent to the spherical harmonics expansion, with the difference 
that the regularity condition can be imposed more easily. 



G(r, 0) = ]T.gn(r) (sin 6) 2n , 



(17) 



H(r, 6) = cos e ( ^2 hn W ( sin ( 



(18) 



The electromagnetic fields are described as 



13, 



Bo 



- sin 9 
arA 



_n=0 

oo 

E 

n=0 



hi 



Inzu 



(sin6») 2 



2 / 

zur g„ 



E r = —r~ 
A 



go + 



71=1 



+ 2(n + l)h n+1 - (2n + l)h n } (s 
(2n- l)h n -i)\ (sin6 



sin cos 
arA 



\g'n + —o(2nh 
^ {zur 2 h' n - 2(n + l)ffn+i } (sin 0) 



where a prime ' denotes a radial derivative and 



= 1 



/ zur \ ' 
V a J 



(19) 



(20) 



(21) 



(22) 



(23) 



Note also that the expansions 1171 and 1181 are not necessarily valid in whole space. The denominators of eqs. 11911 - 1221 
become zero at the light cylinder, which is located at b(r) sinS = 1. Taking the limit of sinS — > l/b(r) in the numerators, we 
have two regularity conditions among the radial functions, g n and h n 



E 

and 



f , t zu 

i 9n + — 



2n - 



2n 



b 2 



= 0, 



n=0 



(h'n 



2nzu 



0. 



(24) 



(25) 



These relations should be satisfied for the region 1/b ^ 1, which is equivalent to r ^ r c , where the radius r c is defined by 
b(r c ) = 1 and, of course, r c is reduced to 1/fi in the flat spacetime. Note that it is not necessary for the radial functions 
to satisfy the conditions 12 It . 12 "it in the inner region r < r c . It may be possible to construct a global solution by smoothly 
matching the radial functions across the radius r c , which are solved with and without the constraints I24H . 125H on either side 
of r c . We however look for the solutions that satisfy (at least approximately) the constraints l|24|l . i12.~)l in the whole space 
exterior to the stellar surface. This method may restrict the class of solutions, but avoids to locate and handle a matching at 
r c . If no solution can be found satisfying the constraints 12 It . I|25|l in in entire region outside the star, we will then look for 
solutions under weaker constraints. 

We consider the electromagnetic field described or approximated by a finite number of functions g n (r) and h n (r). For 
example, the field by rotating magnetic dipole in a flat spacetime can be described by g n , h n (n = 0, 1). See Appendix for the 
results in the flat case for comparison. It is natural to approximate the electromagnetic fields in a similar fashion also in the 
case of a slowly rotating spacetime. 



The function H is not a scalar one, but rather a " pseudo-scalar" . 
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The regularity conditions for the non-vanishing components g n , h n (n — 0, 1) are given by 



= X = g[ + b 2 g' - ^h ? (3 - 2b 2 ) hi 

Or vor^ x ' 

and 

= y = h' 1 + b 2 h -^-g 1 . 

Using these conditions in eas. ll9t - 122L we have 

B r — —h'o cos 6, 

(ho — 2hi — ujr g'o) sind, 

3 / 2 \ 2 
E r = -g'o sin 2 = - (g'o H -hi ) H ^/iiP2(c 

1 2 , 
i?e = — (zjr h — 2gi) sin 8 cos 6*. 



1 

ar 



(26) 

(27) 

(28) 
(29) 

(30) 

(31) 



The electric field is given by a mixture of monopolar and quadrupolar terms. In particular, when the total electric charge is 
zero, the monopolar part vanishes, with the quadrupolar one being the one induced by rotation. This condition can be written 
as 

2 



go - 



rhi =0. 



The divergence of the magnetic field yields 
V--B = -aHcostf, 
where 



H = h'o + -h'o + ^-g'o + z^Z2( 2hl ~ h °> 



(32) 



(33) 



(34) 



Using the condition 1)32(1 in ea. l|34|l . we have a second order differential equation of ho only. Enforcing that the magnetic field 
is divergence-free is equivalent to setting Ti — so that equation 134H has solution 

2M 2 



3/i 

4M 3 



ra 2 In a 2 



+ 2M 



(35) 



where n denotes the magnetic dipole moment. This relativistic solution for the magnetic field is well known in the literature 
jGinzburg fc Ozernov lllflfid : lAnderson fc Cohen lll97fl IPetterson lll97^ . 
We now express the divergence of the electric field as 



V -E = -aQ+Qsm 
where 



r a 2 i 

y = go + -go 

r 



and 



Q 



3o 



2ro , 4 



-hx-y 



(36) 
(37) 

(38) 



In vacuum and in the absence of a net electrical charge, the solution to Maxwell equations is determined by the following set 
of differential equations 



H = Q = , 



(39) 



together with the regularity conditions X = y — 0. It is apparent that this is an overdetermined system if zu'hi/zu ^ 0, for 
in this case the two equations Q = and y — cannot be satisfied in the presence of spacetime dragging. 

An additional constraint hi — leads to a trivial solution, g n = h n = and it is therefore impossible to impose hi = 
everywhere. As a result, the solutions found with this method are only approximate solutions to the Maxwell equations in 
vacuum, with the deviation being measured through the term w'hi/zu. If the value of this quantity, which clearly depends on 
position, is small enough, then the level of approximation can be very good. Hereafter we will monitor the term ■uj'hi/m as a 
consistency check to the solutions found. 
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The functions gi and hi are obtained from Q = X = 0. 

(40) 



gi = ^h' + ^- ci |3r 3 (r 2 - 3Mr + 2M 2 ) In a 2 + 2Mr 2 (3r 2 - 6Mr + M 2 )} + 4j 



hi = ~& (41) 



4r 



3 7 

ci {6r 3 (2r - 3M)lna 2 + 4Mr(6r 2 - 3Mr - M 2 )} + — (rlna 2 + 2M) 



(42) 



M 3 

where a constant c\ is given by the values at surface (r = R) such as q_r = a(R), 

1 (i? 2 In q| + 2A/J? + 2M 2 )R 2 tt - 2J(R In q| + 2M) 
Cl ~~4.RM 3 3R 2 {R- M)a 2 R lna 2 R + 6MR 2 -12M 2 R + 2M 3 ' ( ^ 

We have assumed that the stellar surface behaves as a perfect conductor so that the following boundary condition can be 
imp osed at the stellar surface : Eft + B r a~ zuR sin 6 = 0. Note that these analytic solutions have already been derived elsewhere 
(e.g. iKonno fc Koiimal (2000); R ezzolla et all l)2001alb|^ . 

By direct calculations, we can check that these solutions satisfy Q = 0, but never y = 0. Indeed, these solutions are such 
that y = vj'hi/va ^ and for r S> M this term can be easily estimated to be 

-/ 1 i~-^(4A/ 6 Ci +5J)-^ ? ocr" 6 . (44) 
vj 10 v ' Mr 

Since this term rapidly decreases with distance, the approximate solutions will be rather accurate at a sufficient distance from 
the stellar surface and could therefore be used in practical applications. 

The analytic expressions 135L 140H . 141H and 1421 are the same ones as obtained in previous works, although with some 
slight differences in notations. Note that the "error term" 144t is effectively a second-order term in Q. and is therefore absent 
in previous analyses which were restricted to the low-frequency approximation. Because in our treatment such a term is not 
neglected a priori, we can explicitly check the accuracy of the expressions l|35^ . I|40|l . I|41|l and l|42|l at all radial positions. 
Indeed, these solutions represent a rather good approximation and can be used in practice in a slowly rotating spacetime, if 
the "error term" (which is usually quite small) is within a tolerable range. 



4 NON-AXIALLY SYMMETRIC CASE 

In the previous Section, the axis of the magnetic dipole is aligned with the rotational axis. Here, instead, we consider what 
happens when the dipole is misaligned. In particular, we focus on the case in which the axis of the magnetic dipole is 
perpendicular to the rotational axis. 

The method used here to solve the Maxwell equations is almost identical to the one used in the previous Section for the 
axially symmetric configuration. In order to point out the similarities, we will adopt the same notation for the radial functions 
g n ,h„. While this may induce some confusion, it is sufficient to bear in mind that the basis functions used in the expansion 
are different in the two cases [cf. eqs. (117P — J18| and 1451 - 146^ , 

Since the potentials G and H would depend on the azimuthal wave number m — 1, they are expanded by the spherical 
harmonics Y 2 * an d Yn+ii respectively. As a result, the angular expansion can be written as 

G(r,6,<t>,t) = e i<6 ' sin (9 cos (9 ^^ ffn (r)(sin(9) 2n ^ , (45) 



H{r,6,<j>,t) =e i</, 'sin0 |^ h„{r) (sin 6) 2n J , (46) 
where 

4>' = 4>- fit. (47) 

The functions should depend only on r,9,(f>' from the stationary condition as mentioned in Section 2. As shown in eas. lim 
and 1111 . the denominators in the poloidal electric and magnetic fields become zero at the light cylinder. Using the expansion 
forms, we have two regularity conditions, which should be satisfied for r ^ r c . The conditions for the radial functions are 
written as 



E 

n = 



{g' n + (2n+l)^h n }b- 2n =0, (48) 
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EK-5( 2n+1 -ir)+" 2 "= 



(49) 



In a spherically symmetric and static spacetime the electromagnetic fields are described completely in terms of four radial 
functions g„,h n (n = 0, 1). This set of functions can also provide approximate expressions in the case of a slowly rotating 
spacetime. The regularity conditions for the non- vanishing components are given by 



zub 2 



3zu , 



= X = g[ + b 2 g' + —ho + 

a 2 ur 

o = ys / l ; + 6 X-5(6 2 -%o-5(3-l) ff , 

Using these conditions, the expressions for the electromagnetic fields regular at the light cylinder are 

B r = 



, w 

"o 2 9o 

or rur 



4 • 2, 

2 5i sin I 



sinf 



Be = — — \ho + (b 2 h() + 3hi + zur 2 g' ) sin 2 0] cos 0, 

OLT J 

B<t, = — — [ho + hi sin 2 6\ e l< * , 



9o + -^h 



E 6 = -— [go - {(2 - b 2 )g a - 3gi + zur 2 h' } sin 2 9] e^' , 



E,-> 



ar 

i 



[go + gi sin 2 6\ cos ( 



The divergence of the electric and magnetic fields yields 
V • E — —aQ sin 6 cos 9, 
V-B = -aHsm9- Qsin 3 6», 
where 

r ct (ar) 2 (ar) 2 r 2 \ a 2 " 



H 



Q = 



fto + -h + —go - — I — — I 9o + 

4a ( zu 1 



cr r 2 \ a 2 ) (ar) 2 



h^ + ^-^Cib 2 -2)h , 



(ar) 2 



zur^ \ zu 



-gi - x 



(50) 
(51) 

(52) 
(53) 
(54) 
(55) 
(56) 
(57) 

(58) 
(59) 

(60) 
(61) 
(62) 



The solutions g n ,h n (n — 0, 1) of the Maxwell equations in vacuum are determined by Q = TL = Q = Q and the regularity 
conditions X = y = 0, but it is clear that the system of differential equations is overdetermined in a slowly rotating spacetime 
because of the presence of a redundant equation. In this respect, and as mentioned in the previous Section, the solutions to 
these equations are not exact. 

In order to obtain the Deutsch-type solution, we impose the following conditions (see the Appendix for details) 



gi = 0, hi = -- (zur 2 g' + b 2 h ) 



(63) 



Using these conditions, we have X — and Q — 0. The two functions go, ho are determined by solving Q = H = 0. We 
introduce the functions Fi(l = 1, 2) so as to express ho and go as 

5() = 1T F2 ~! Fl ' (64) 



a . zu 
= ~2~ + IT ' 



where .Fi and i<2 satisfy the coupled equations 
a 2 (a 2 Fl)' + (zu 2 -^\ Fl + \a 2 zu'F 2 = 0, 



<<-[<<- { ) :!n 2 -'.F, = it. 



(65) 

(66) 
(67) 
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It is easy to check that the functions 1641 and 1651 satisfy the coupled equations Q = Ti — with the conditions 1631 . Indeed, in 
the limit of a non-rotating sp acetime, equations (1661 and 1671 1 become the Regge- Wheeler wave-like equations with frequency 
Q for I = 1,2 JwheeleJll955n . If the spacetime is rotating, however, the two equations are coupled because of the spacetime 
dragging through the coefficient za' = -co' — bJ/r 4 . From this choice of functions, we have 

y = ~F 2 . (68) 

This does not vanish in general, but rapidly decreases to zero with the radius because of the factor 2J/r 4 . In fact, we can find 
y ~ -6nJQR 2 /r 6 near the light cylinder and y ~ 2/iJQ 3 R 2 exp{iQ(r + 2M\n{r/2M - l))}/r 4 beyond the light cylinder. 

Since the function F2 describes the induced electromagnetic fields by a rotating perfect conductor, it is a first-order 
quantity in Q. Consequently, the right-hand-side of 16811 is second-order in Q and was thus absent in previous works performed 
in the slow-frequency approximation. Furthermore, eqs. 1661 and 1671 are sig nificantly simpl i fied wh en the terms 0(fl 2 ) can 
be neglected and, in fact, analytic expressions have been found in this case bv lRezzolla et al l i2001allbh . It should be stressed, 
however, that such analytic solutions can be valid only in the vacuum exterior close to the stellar surface, i.e. for r < r c . 
Beyond the light cylinder, in fact, the wave-like nature of the electromagnetic fields becomes important and the term 0(Q 2 ) 
in eqs. 1661 and (16711 can no longer be neglected. 

Expressing the electromagnetic fields in terms of functions Fi we have 

Vi sin ee»' , - f £ F{ + -^F 2 ) cos 9e* , - i ( £-F[ + ^F* cos 26) 
V s \2r bar ) \2r bar ) 



B = 
E = 



e 



-J_f 2sin2 (9e^', f^-Ft-— F 2 'cos26)e l4 '\ i {—Fx - — Fi] cosfle 1 ^' 
2r A \ 2ar or / V 2ar br ) 



(69) 
(70) 



and we can use the perfect conductor condition at the surface, i.e. Eg+B r a 1 vjR sin 6 — 0, which is reduced to F£ —3ra7.Fi /a 2 = 
at r = R. 

Using this boundary condition, the numerical solution of the coupled equations 1661 and 1671 for the magnetic field B r 
is shown in Fig.l for a typical choice o f the st ellar parameters. In the same figure, we also show the approximate solution 
in cur ved spacetime by iRezzolla et al ] (hOQlallrl [indicated by (a)] and the asymptotic one in flat spacetime by iDeutschl 
lll955l) [indicated by (f)]. The function B r decreases with radius and the difference between our numerical result and analytic 
expressions is not so clear when shown in absolute magnitude. To highlight the differences we have therefore multiplied the 
magnet ic field intensity B r in Fig.l by the factor r 3 /(fi sin Qe 1 ^ ). Note that our numerical solution agrees rather well with 
that of lRezzolla et al ] feooiall t ^ near the sur face. Around r ~ 0.2S1 , our result deviates from it, approaching however the 
solution for flat spacetime by Deutscb| (Il955l) as the distance from the stellar surface increases. 

The explanation for this behavior is rather simple since the analytic expression for curved spacetime in the low-frequency 
approximation is valid only well within the light cylinder region r ~ 0.2O" 1 < r c , where the relativistic corrections are large 
and the electromagnetic frequencies small. Indeed, in the low- frequency approximation the point r c ~ 1/fi is regarded as 
"infinity" and, as expected, the results of the low-frequency approximation c annot be accura t e near this limit. As the distance 
from the star increases, a qualitative difference appears from the solutions of lRezzolla et al 1 (bOQlallrl and this is represented 
by the sinusoidal, wave-like nature which is clear beyond the light cylinder. In this region, the gravitational corrections are 
not important and, in fact, our numerical result agrees well with that of flat spacetime. Fig.l therefore serves to illustrate 
how our approach is very useful in connecting two different regimes, i.e. the relativistic one and the and the asymptotically 
flat one within a single, semi-analytic approach. 



5 CONCLUDING REMARKS 

We have studied the stationary solutions of the Maxwell equations in slowly rotating spacetime. The electromagnetic fields 
exterior to rotating magnetic dipole moment have been derived with the perfect conductor condition at the surface. This 
problem is not new, but it has been here considered through a different type of approach and with a different level of 
accuracy. In particular, we have considered the general relativistic electromagnetic fields exterior to a rotating neutron star in 
the slow-rotation approximation but not in the low-frequency approximation. The solutions derived here reduce to the exact, 
analytic ones if the spacetime is spherically symmetric and no frame dragging effect is presen t. Furthermore , in this same 
limit, they can be expressed in terms of analytic functions in the limit of flat Minkowski space JPeutschlll955l) . 

Using these results as a guide, we have looked for solutions that have similar angular dependence but more complicated 
radial functional behavior. Adopting a semi-analytical approach in which simple expressions for the electromagnetic fields 
need to be completed through the calculation of coefficients coming from the numerical solution of two partial differential 
equations, we have found approximate solutions to the Maxwell equations. The error terms present in our solutions are O(0. 2 ) 
and were therefore absent in previous approaches using the low-frequency approximat ion. The numerical so lutions found 
in this way agree well with the known analytic results both in the strong-field region jRezzolla et al.ll2001allrl and in the 
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Figure 1. Magnetic field B r r 3 /([i sin Be 1 ^ ) as a function of radius r/r c . The left and right panels show the behavior near the surface 
(R ^ r Sj r c /2) and that for wave region (_R ^ r sC 10r c ). The solid line with symbol 'n' represents our result solved for 2M/r c = 
0.01, R/r c = 0,05, J = 1 0~ 3 r^. The dotted lines wit h symbols 'a' a nd 'f respectively represent the analytic result for curved spacetime 
jRezzolla et al.lboOlallf) and that for flat spacetime fDcutsch 1 9551) . 



asymptotically flat one l|DeutscbJ ll955V where the wave nature of the solution is important beyond the light cylinder in the 
case of non-axially symmetric configurations. Furthermore, our numerical solutions represent a useful tool to determine when 
these two analytic approaches cease to be accurate. Overall, the new method proposed here offers the advantage of being able 
to bridge in a single approximate but accurate formalism both the strongly relativistic and the asymptotically flat regimes. 

Our approach is fundamentally based on a series expansion for which we discuss here only the lowest-order terms. However, 
it is in principle possible to increase the accuracy of our results by increasing the number of expansion series. For example, by 
including the higher-order functions t/2, hi, the solutions can be calculated for a larger system of partial differential equations. 
Also in this case, there would be a redundant equation in the set of equations signalling terms which are incompatible within 
a set of differential equations. While we will not show any explicit calculations here, it is not difficult to realize that in the 
lowest extension to the results presented here the new correction would come through a term ~ vj'K^Ivj (or zu'g2/zu) in the 
system of the functions g n ,h n (n — 0, 1,2). The function hi (or gi) is expected to be small, because it should be zero both 
at the surface and infinity and the overall accuracy is therefore improved in this enlarged system of equations. By increasing 
further the number of functions, the difference between the incompatible equations becomes small. Of course, the limit of 
n^> 1 may be preferable because of its increase in accuracy, but it also has clear limitations in practical applications. 
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APPENDIX A: DEUTSCH SOLUTION 

We here summarize the Deutsch solution and offer it for comparison. The solution describes the electromagnetic fields produced 
by a dipole magnet rotating with angular frequency fi. The quadrupolar electric field is induced by a perfect-conductor 
condition at the surface, written as Eg + QR sin 6B r = at r = R. We here assume that the angle between the dipole moment 
and rotational axis is x so tna t the components of the electromagnetic fields are given by 

B r = ^ cos x cos # + ^jSie !A sin x sin ^, (Al) 
Be = cos x sin # — -^-(s2 — diSs)e lX sinxcos6>, (A2) 
B^ = - ^ (S2 - diss cos 29)e iX sin x, (A3) 
E r = -^^cos X P 2 (e) + ^ S3e ' A sinxsin2e, (A4) 

?2 



Eg = 
= 

where 



2uf27T . „„ /ufi 2nd! A . ix /Ar s 
— A — cosxsm26' - I — si + -^^s 4 cos26< \ sinxe , (A5) 

-i(^si + ^s 4 )sinxcos0e iA , (A6) 
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(A8) 

(A9) 
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